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Abstract. The Radon transform IZf of functions / on SO(3) has recently 
been applied extensively in texture analysis, i.e. the analysis of preferred 
crystallographic orientation. In practice one has to determine the orientation 
probability density function / G L 2 (SO(3)) from TZf g L 2 {S 2 X S 2 ) which is 
known only on a discrete set of points. Since one has only partial information 
about IZf the inversion of the Radon transform becomes an ill-posed inverse 
problem. Motivated by this problem we define a new notion of the Radon 
transform IZf of functions / on general compact Lie groups and introduce 
two approximate inversion algorithms which utilize our previously developed 
generalized variational splines on manifolds. Our new algorithms fit very well 
to the application of Radon transform on SO(3) to texture analysis. 
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1. Introduction 

The objective of the present paper is to introduce Radon transform on compact 
Lie groups and to show how spline interpolation can be used for approximate in- 
version of such transform on general compact Lie groups and in particular on the 
group of rotations SO(3). 

The Radon & transform on a compact Q Lie group associated with a closed 
subgroup assigns to a smooth function / integrals over submanifolds of Q which 
have the form x3ft?y~ x where x, y £ Q. An important fact is that the Radon 
transform S%f is a function on £7 / J^IC x Q / and not merely on Q x Q as one could 
expect (see below for details). The typical problem is to reconstruct / knowing £%f. 
In practice one has only integrals over manifolds {M-^}^ from a finite subfamily. In 
this situation inversion becomes an ill posed problem. Our objective is to consider 
appropriate regularization of this ill posed problem and to develop a method for 
approximate inversion of the Radon transform. 

In our approach to approximate inversion of the Radon transform (section 5) 
we consider inversion as an interpolation problem. Namely, given a set of integrals 
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of an / over a finite family M. = {A4„}^ of submanifolds, we find a "smoothest" 
function which has the same set of integrals as / over submanifolds from the family 
M. 

Our work is inspired by recent applications of the Radon transform on 50 (3) to 
texture analysis, i.e. the analysis of preferred crystallographic orientation. 

The orientation probability density function (ODF) representing the probability 
law of random orientations of crystal grains by volume is a major issue. In x-ray or 
neutron diffraction experiments spherical intensity distributions are measured which 
can be interpreted in terms of spherical probability distributions of distinguished 
crystallographic axes. A first mathematical description of the problem was given 
in 0. 

One equips a crystal with an inner orthogonal coordinate system {ei,e2,es}. 
Additionally one distinguish an outer orthogonal coordinate system {ui, 1x2,143} 
related to the specimen. The orientation of a crystal in the specimen is defined 
by the unique rotation 7 G 50(3) which maps the inner coordinate system to the 
outer one, i.e. gei = m for i = 1,2,3. Note that in this model we neglect the 
spherical symmetries of the crystal. 

The function of interest is the orientation density function ( ODF) f G L 2 (50(3)) 
that is the probability measure on 50(3). Hence the function value /(<?) gives the 
portion of crystals in the specimen with orientation g. 

The practical measurement sends a beam through the specimen coming from 
the direction h £ S 2 and measures the intensity, emitted from the specimen in the 
direction r £ S 2 . One can interpret the result as the integral over all orientations 
g 6 50(3) with g-h = r;h,re S 2 . The set C Kr = {g £ 50(3) : g-h = r;h,re S 2 } 
of those orientations is called a great circle in 50(3). 

Since 50(2) is the stabilizer of £0 G S 2 , where £0 is the north pole one has 

(1) C h , r = h'SO{2)r'- 1 := {h'gr'- 1 , g £ 50(2)} h' , r' G 50(3), 

where h',r' G 50(3) satisfy h' ■ £0 = h and r' ■ £0 = r - Here (x,y) is a fixed point 
in 5 2 x S 2 . 

Definition 1.1. The Radon transform of a continuous complex-valued function f 
on 50(3) is a function on S 2 x 5 2 which is defined by the formula 

(2) af(x,v)= [ f(g)dx. 

This transform can be extended to all functions in L 2 (50(3)). 

From the very definition Z%f is a function on 5 2 x 5 2 . Moreover, it is shown 
that £%f is in the kernel of the Darboux-type differential operator (also called ultra- 
hyperbolic operator). In other words, one has 

A x <%f(x, y) - A y ^f(x, y), (x, y) G 5 2 x 5 2 . 

It is interesting to note that this condition is similar to the classical condition of F. 
John I5j for A-ray transform in M 3 . 

The Darboux-type equation shows that fflf belongs to the span of the tensorial 
product of spherical harmonics y].yi which is obviously a subspace of L2(S 2 x 5 2 ) 
whose basis is y l k yf ■ 

Because 5 3 is a double covering of 50(3) there is also a close connection to the 
spherical Radon transform of even functions considered by S. Hclgason ([18], [19] ) 
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as a special case of Radon tranforms on homogeneous spaces. A comparison of 1% 
to the spherical Radon transform of even functions is given in [4]. 

The numerical aspect of solving the crystallographic problem as well as the 
inversion of is treated in [6] and [20] . 

To tackle the real-world problem of sharp textures, i.e. the pole figures consists 
mainly of a few delta peaks, wavelets on the sphere S 3 and the group SO(3) and 
their behavior under had been studied in [3] and [10]. In [TT] and [10] an 
grouptheoretical approach is used to generalize the results for wavelets on groups 
and homogeneous spaces which leads to a new class of wavelets, called diffusive 
wavelets. 

Another approach to solve the inversion problem of is done in [8] , where Gabor 
frames on the Spin(3) had been used to solve the problem numerically. 

In section 2 of the paper we use a group theoretic and representation theoretic 
approach to describe the Radon transform on a general compact Lie group Q. At 
first hand f){x,y) seems to be defined over Q x Q while deeper inverstigations 
reveal that is invariant under right shifts of x as well as y and hence 3$ is 
rather defined over Q / J^f x Q/Jtf'. For the practical application Q = SO(3) and 
,ye = SO(2) and thus Q/J^ = SO(Z)JSO{2) = S 2 . In this particular situation a 
number of explicit inversion formulas is known. For more details we refer to [2] , [I] , 

0, H, [H], m, |2a- 122| . 

In section 4 of the paper we discuss generalized variational splines on compact 
Riemannian manifolds (see [30], [34], [33]) and apply this concept to the Radon 
transform. 

In [30]- [34] a theory of generalized (=average) interpolating variational splines 
was developed in compact and non-compact manifolds of bounded geometry. We 
use term generalized splines to stress that we interpolate functions on manifolds by 
using values of their integrals over submanifolds of a given manifold. It includes the 
classical way of interpolating functions is by using their values on discrete sets of 
points. Variational splines on manifolds were used in our papers for reconstruction 
of Paley- Wiener functions on manifolds. In [31] the theory of generalized varia- 
tional splines was applied to the spherical Radon transform of even functions (Funk 
transform) and to the hemispherical transform of odd functions on n-dimensional 
spheres. 

The development in the present paper is slightly different from our approach in 
previous papers in the sense that here we introduce splines using general elliptic 
second order self-adjoint operators and not only the Laplace-Beltrami operator. 
We prove existence and uniqueness of interpolating variational splines associated 
with a general elliptic second order self-adjoint operator on a compact Riemannian 
manifold. We remind (see Theorem I4.6j) that every generalized spline is a lin- 
ear combinations of what we call generalized fundamental solutions (= generalized 
Green functions). Using these results "explicit" formulas of variational splines in 
terms of eigenfunctions of the corresponding operator are given. 

In section 3 we develop our first method for approximate inversion of Radon 
transform on compact Lie groups. Note that in a similar situation for Radon a 
hemispherical transforms on unit spheres this idea was introduced in [34] . The 
idea is that one can always have an approximate inversion of " any kind" of Radon 
transform by treating inversion as a generalized interpolation problem. Thus, we 
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invert group Radon transform by "interpolating" a function of the group using its 
integrals over a finite family of appropriate submanifolds. 

In the last section we review our fundamental inequalities (Lemma 16.21 Theorem 
16. 3[) for functions with many zeroes (see (30) . |34) . [33]) and our Approximation 
Theorem for interpolating variational splines on general compact manifolds ([34], 
[33]). The rate of approximation in this theorem is expressed in terms of Sobolev 
norms. It is worth to note that using standard methods of interpolation theory and 
the fact that Besov spaces on manifolds are interpolation spaces (real interpolation 
method) between two Sobolev spaces the corresponding inequalities can be extended 
to Besov norms. 

We use these results to describe our second method of approximate inversion of 
the Radon transform on SO(3). The idea of this method is to interpolate not the 
function but its Radon transform on the manifold S 2 x S 2 and then "to return" 
to the original manifold SO(3). This method allows to obtain explicit estimates of 
degree of approximation in Theorem 6.5. Moreover, in the same section we establish 
a Sampling Theorem for Radon transform on 5*0(3) (Theorem 6.6) which says that 
one can have a complete reconstruction of w-bandlimitcd (on SO (3)) function / 
by using only the values of its integrals over sufficiently dense specific family of 
submanifolds of SO (3). 

It should be stressed that the variational spline problem fits in some sense optimal 
to the practical question of determining the ODF / from measurements of the 
Radon transformed /. On the one hand we are interested in regions where the values 
of / are large but if the curvature of / is small in those regions it would be more 
useful to increase the measurements around the maximum of / and those points 
where the curvature is large. Hence the right criteria to increase the density of 
measurements around points where (1 — A)/ is large. The density of measurements 
should be high at those points where the interpolation is highly nonlinear. 

It should be noted that a very different and more constructive approach to splines 
on two-dimensional surfaces was developed in pQ, [9], [12]- [14], [23] . 

An approach to splines on compact manifolds which is similar to our approach in 
[30]- [34] (in the sense it is based on a "zeroes lemma" and Green's function repre- 
sentations) was recently developed in [16] , [17] for the classical way of interpolation 
on discrete sets of points. 



2. Fourier analysis and Radon transform on compact Lie groups 

2.1. Fourier Analysis on compact groups. One of the most important theo- 
rems of functional analysis is the spectral theorem for compact self-adjoint operators 
on a Hilbert space. Which states that if A is a compact self-adjoint operator on a 
Hilbert space V, then there is an orthonormal basis of V consisting of eigenvectors 
of A and each eigenvalue is real. The analog of this theorem is the Peter- Weyl 
theorem. We want to recollect some basic notations and properties. 

Let Q be a compact Lie group. A unitary representation of Q is a continuous 
group homomorphism n: Q ^ U{<1t T ) of Q into the goup of unitary matrices of 
a certain dimension d v which will be explained later in the Peter- Weyl theorem. 
Such a representation is irreducible if ir(g)M = Mir(g) for all g G Q and some 
M G <£ d T y - d T implies M = cl is a multiple of the identity. Equivalcntly, C^" does 
not have non-trivial 7r-invariant subspaccs V C C " with n(g)V C V for all g G Q. 
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Two representations m and 7T2 are equivalent, if there exists an invertible matrix 
M such that iri{g)M = M-K 2 {g) for all g eQ. 

Let Q denote the set of all equivalence classes of irreducible representations. 
Then this set parametrcrizes an orthogonal decomposition of L 2 (Q). 

Theorem 2.1 (Peter- Wcyl, 35 ). Let Q be a compact Lie group. Then the following 
statements are true. 

a: Denote H n = {g >->• trace(7r(p)M) : M G C d " xd "}. Then the Hilbert space 
L 2 {Q) decomposes into the orthogonal direct sum 

(3) L 2 (C7) = 0H. 

b: For each irreducible representation 7r G G the orthogonal projection 
L 2 {Q) — > %-n is given by 

(4) f^d„ f{h) X 7r(h' 1 g)dh = d n f * X ^, 

Jg 

in terms of the character x^(g) = trace(7r(g)) of the representation and dh 
is the normalized Haar measure. 

We will denote the matrix M in the equation / * \-k = trace(7r(g)M) as Fourier 
coefficient f(n) of / at the irreducible representation tt. The Fourier coefficient can 
be calculated as 

/>)= / f(g)n*(g)dg. 

Jg 

The inversion formula (the Fourier expansion) is then given by 

f(g) = ^trace(7r( 5 )/(7r)). 

ttEG 

If we denote by ||M||^- S = trace(M*M) the Frobenius or Hilbert-Schmidt norm of 
a matrix M, then the following Parseval identity is true. 

Corollary 2.1 (Parseval identity). Let f G L 2 (Q). Then the matrix-valued Fourier 
coefficients f G <C d7,xd7 ' satisfy 

(5) \\f\\ 2 = J2 d *WfW\\HS- 

tt£G 

On the group Q one defines the convolution of two integrablc functions /, r G 
L\g) as 

f*r(g) = [ f(h)r(h- 1 g)dh. 
Jg 

Since / * r € L 1 (Q), the Fourier coefficients are well-defined and they satisfy 

Corollary 2.2 (Convolution theorem on Q). Let f, r G L 1 (Q) then f * r G L 1 {Q) 
and 

f * r(?r) = /(7r)f(7r). 

The group structure gives rise to left and right translations T g f i— > fig -1 -) and 
T 9 f >-> f{-g) of functions on the group. A simple computation shows 

T)J{n) = fay*(jg) and fZ~f(ir) = n(g)fa). 



(i 



SPLINES FOR RADON TRANSFORM ON COMPACT LIE GROUPS 



They are direct consequences of the definition of the Fourier transform. 

The Laplace-Beltrami operator Ag on the group Q is bi-invariant, i.e. commutes 
with all T g and T 9 . Therefore, all its eigcnspaces are bi-invariant subspaces of L 2 (Q). 
As "Hjr are minimal bi-invariant subspaces, each of them has to be eigenspace of Ag 
with corresponding eigenvalue — Aj. Hence, we obtain 

Agf = - E A - trace(7r(<?)/(7r)). 

Motivated by situation which was described we introduce a new type of Radon 
transform which is an abstract version of the crystallographic Radon transform. 

2.2. Radon transform on compact groups. Now, we are able to define the 
Radon transform. 

Definition 2.2. Let Jif be a closed subgroup of the compact Lie group Q . The 
Radon transform of a continuous function f £ C(Q) is defined by 

(6) &f(x,y)= [ fixhy-^dh x,yeG, 

where dh here is the normalized Haar measure on Jif. 

Next, we explain that M maps Q into Q / Jtf? x Q / 3^ . For that we use the averaging 
method. 

For the following discussion we mention, that functions on Q / Jf? can be regarded 
as functions on Q, which are constant over right co-sets of the form gjtf? = {gh, h £ 
J^} for all g S Q. The projection of a function on Q onto functions on Q j ffl 
corresponds to an averaging method over gffl: 

(7) Pjr/Cff) = / f(9h) dh. 

It can be shown, that Wjp in Fourier domain acts by multiplying the Fourier coef- 
ficients f(ir) by 

(8) tt^= / n{h)dh 

from the right. Further njg> is a projection and without loss of generality ttj? = 
diag(l, 1, 0, .., 0), where the number of l's corresponds to the number of 
invariant vectors in the representation Hilbert space of w. For details on the pro- 
jection see [3B]. Next we discuss the Range of ffl. Since x, y in ^ are elements of 
Q at a first look it seems that the Radon transform is defined over Q x Q. While 
a deeper investigation reveals that &f(x,y) is invariant under right shifts of x as 
well as under right shifts of y, hence Si is rather defined over Q j ■Jtf' x Q/J^f. 

Lemma 2.3. The Radon transform M maps functions over Q to functions over 

g/,ye x g/je. 

Proof. We look at S% on the Fourier domain. Let first y G Q be fixed and regard 
^fi'iV) t° De a function on Q in the first argument, then 

(9) @K^y){v) = -K X K*{y)f{-K)-KeQ, 

Hence the the function Mf(-,y) is invariant under the projection ¥jg>, since the 
Fourier coefficients are invariant under the left multiplication by 7r#»: 

7r^7r^7r*(y)/(7r) = 7r^7r*(y)/(7r). 
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Consequently, 

(10) &f(x-h,y)=&f(x,y) hEJf. 

A look at the Radon transform as function in the second argument y, while the 
first one x = Xq is fixed, we find 



F^Mf(x ,y)= / &f(x ,yh)dxh = 
Jje 

/ d 7r (f(Tr)Tr(x ))Tr.jtfiT(h~ 1 y~ 1 ) dxh 

Tree 

= ^2 d^(f(Tr)TT(x ))7TjifTT*(y) = &f(x ,y), 



Hence Mf(x, y) is constant over fibers of the form yjf? also in the second argument 
and 



(11) ^/(z,0(7r)=7r^7r*(x)/(7r)*, 

where the complex conjugate of the matrices is taken componentwise. Consequently 
^ maps functions over Q to functions over Q / 'Jff x Q/Jtf. □ 

Below we discuss the Radon transform of the space L 2 (Q), its range and its 
inversion. 

Lemma 2.4. Let J$? be the subgroup of Q, determining the Radon transform on Q 
and let Q\ C Q be the set of irreducible representations with respect to Jf. Then 
for f e C°°(g) it is 



(12) \\&f\\l>w/jexo/je) = E rank(7r^)\\m 



\hs- 



Proof. For the proof we expand &f(x,y) for fixed y as function in x over Q (or 
better Q /Jti?) and apply Parseval's identity ([S]), with © we have 

\\^f\\h{g/jexg/M>) = Yl d ^ i IK^^/Mllffs d V 

J Q 

wee 

= X] dw y tracc (/*( 7r ) 7r (y) 7r -#' 7r *(y)/( 7r )) d y 

Tree e 

= ^ d w trace ^/*(tt) ^ Tr(y)nj%>Tr* (y) dy f(ir) 



rree 

= ^ rank(7T^.)trace(/*/) = ^ rank(/T^>)||/(7r) 
•rreei ireei 

Where we made use of 

/ rankTr^f 

y Tr(y)irjen*(y) dy = I ^ y nik(y)n k j(y) dy 



\hs- 



k=l 



(13) = rank(7r ^ } M 

flrr 



□ 
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3. Harmonic analysis and Radon transform on SO (3) 

3.1. Special functions of SO(3). For the Hilbcrt space L 2 (S 2 ) we use the or- 
thonormal system of spherical harmonics in the {yj., i = 1, ...,2k + l,k G No}. 

In case of SO (3) the situation is very comfortable, since every irreducible rep- 
resentation is unitary equivalent to a irreducible component of the quasi regular 
representation in L 2 (S 2 ), given by 

(14) T(g) : /(£) >-» f{g~ x ■ x), 

where • denotes the canonical action of SO (3) on S 2 . The irreducible invariant 
components of L 2 {S 2 ) under T are T-Lj. = {3^, i = 1, 2fc+l}- spanned by spherical 
harmonics of degree k. T k shall denote the irreducible representation, obtained by 
restriction of T to Hk- The matrix coefficients of T k are the Wigner polynomials 
T*. of degree k: 

2fe+l 

(is) yiig- 1 ■ o = J2 T !sM(Z) T k ( 9 ) = {yi{g- l -),yi{-)) LHS *y 

By construction the projection matrix on the Fourier handside given in (|5]) is 
7r SO(2)(^) = diag(l, 0, 0) where number of zeros is 2k and n SO (2) (k) of dimension 
(2k + 1) x (2k + 1). Since matrix coefficients always have the norm 4-, where is 
the dimension of the representation, we have 

(16) T*(g) = Jj^y k (g.S ), 

where £o € S 2 is the base point of SO(3) / SO(2) ~ S 2 , often chosen as north pole 
and its stabilizer is the factorized subgroup SO(2). 

The eigenvalue of Laplacian on 5*0(3) and on S 2 corresponding to polynomials 
of degree k is — k(k + 1): 

(17) A so{3) T k = -k(k + i)T k , A S2 yi = -k(k + i)yl 

Definition 3.1. The subgroup Jtf is called massive, ifr&iikJt?(ir) < 1 for all ir G Q. 

Furthermore, an irreducible representation tt G Q is called class-1 representation 
with respect to the subgroup , if rank,34? (tt) > 1. 

Lemma 3.2. ( 36, Chapter IX. 2. 6]) SO(n) is a massive subgroup of SO(n + 1). 
Furthermore, the family T k , k G No, gives up to equivalence all class-1 representa- 
tions of SO(n + 1) with respect to SO(n). 

For the following we fix the 'north pole' £o of S 2 . Then the set of zonal spherical 
harmonics is one-dimensional and spanned by the Gegenbauer polynomials. Further 
the dimension of zonal functions in Hk is one for all k > and it is spanned by 

Gegenbauer polynomial of order C fc 2 (£o • £), note that cos(Z(£o, £) = Co • 0- Hence, 
to be zonal on S 2 for a function /(£) means to be invariant under the action of 
SO(2), i.e. to depend only on the angle between the argument £ and the 15*0(2) 
invariant point £o- 

The following addition theorem holds true 

Theorem 3.3 (Addition theorem). For all £,77 G § 2 and k G No 
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3.2. Radon transform on 50(3). The Radon transform on 50(3) is defined in 
Dcfimtion ll.il We also need apropriate function spaces, which adjust the smooth- 
ness of the functions. Below we define Sobolev spaces for S 2 x S 2 . 

Definition 3.4. The Sobolev space H t (S 2 x 5 2 ), t 6 R, is defined as the domain 

t 

of the operator (1 — 2As2 x s^) 2 with graph norm 

\\f\\ t = ||(1 - 2A 52xS2 )§/|| i2(s2xS2) , / e L 2 (S 2 x 5 2 ). 

Since 

(19) @{T k ){x,y) = T k {x)n so{2) {T k {y))* 

we have 



(20) &T%(Z,ri) = T&WTfcr,) = ^-^(^fo). 

This formula shows that range of Si belongs to kernel of the Darboux-type operator 
i.e. 

(21) A x Mf(x,y) = A y <%f(x,y), f e L a (SO(3)). 

Because fflf is in the kernel of the Darboux-type operator we also need the following 
Sobolev-type function space. 

Definition 3.5. The Sobolev space £T t A (5 2 x S 2 ), t € R, is defined as the subspace 
of all functions f 6 H t (S 2 x S 2 ) such Aif = A 2 f. 

Now we define Sobolev spaces on SO (3). 

Definition 3.6. The Sobolev space H t (SO(3)), t G R, is defined as the domain of 

t 

the operator (1 — 4Aso(3)) 2 with graph norm 

HI/HI; = ||(1 - 4A so(3) )3/|| i2(so(3)) , / e L 2 (SO(3)). 

Because the operators 1 — 2A s2xS2 and 1 — 4A5 ( 3 ) have positive spectrum 
and one can choose corresponding eigenfunctions which form an orthonormal basis. 
These definitions are consistent with another norm on the Sobolev space H t (M) 
which will be given in Definition 14. II 

Theorem 3.7. (Range description) For any t > the Radon transform on SO(3) 
is an invertible mapping 

(22) ^ : H t (SO(3)) ^ HfiJS 2 x S 2 ). 

I+ 2 

Proof. It is sufficient to consider case t = 0. If = 2k + 1 is the dimension 
of the irreducible representations and —X 2 = —k(k + 1) are the eigenvalues of 
the Laplacian Aso{3) we have dk = \Jl+ 4A 2 ,. Since J^f = SO{2) is massive in 
Q = 50(3) and T k are class- 1 representations of 50(3) with respect to 50(2) we 
have rank(7Tj?r) = 1 and Gi = G- Now the assertion follows from (fT2]) and ([T3|) . Wc 
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have 



Il^/Hl = ||(1 - 2A s2xS2 )i^/||| 2(s2xS2) = Y^d* [ \\Vd^^*(y)f^)\\ 2 HS dy 
2 Jg 

neg 

= Y d l J tracc (/*( 7r ) 7r (y) 7r ^ 7r *(y)/( 7r )) d 2/ 

Tree e 

= E 4 tracc y 7r(y)7r.^7r*(y) dy/(7r 



wee 

= ^ ^rank(7r^)trace(r/) = ^ drll/WHlrs = I l/l \h(SO(3)) = 1 1 1 1 1 

□ 

From Theorem 13.71 we deduce the reconstruction formula for the Radon trans- 
form on 5*0(3). This is also the result of a simple calculation involving spherical 
harminoncs as well as Wigner polynomials. Using the identities (|20[) we get the 
following theorem: 

Theorem 3.8. (Reconstruction formula) 
Let 

oo 2fe+l 



(23) (®g){x,y) = f(x,y) = ]T £ f{k,i,j)y k {x)y> k (y) € Hf(S 2 x S 2 ) 

k=0i,j=l 

be the result of a Radon transform. Then the pre-image g G L 2 (SO(3)) is given by 

oo 2/c+l /„, i 1 \ 00 

( 24 ) # = E E { -^ 2 f(k^,j)T^=Y J ( 2k + 1 )^(9(k)T k ) 

k=0 i,j=l k=0 

(25) g(k,i,j) = ^-f(k,i,j). 

4. Generalized variational splines on compact Riemannian manifolds 

We consider a compact Riemannian manifold M without boundary. Let £ be a 
differential of order two elliptic operator which is self-adjoint and negatively semi- 
definite in the space L2(M) constructed using a Riemannian density dx. The 
spectrum of such operator always contains Ao = 0. In order to have an invertible 
operator we will work with I — C, where / is the identity operator in L%(M). It is 
known that for every such operator C the domain of the power (— Cfl 2 , t G K, is 
the Sobolev space H t (M). There are different ways to introduce norm in Sobolev 
spaces. We choose the following definition. 

Definition 4.1. The Sobolev space H t (M),t G IK can be introduced as the domain 
of the operator (1 — C) 1 ! 2 with the graph norm 

\\f\\ t =\\(l-Cy/ 2 f\\,f€H t (M). 

Note, that such norm depends on C. However, for every two differential of order 
two elliptic operators such norms are equivalent for each tsl. 

Since the operator (— C) is self-adjoint and positive semi-definite it has a discrete 
spectrum = Aq < Ai < A2 < and one can choose corresponding eigen functions 
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ipo,(pi, ... which form an orthonormal basis of L2(M). A distribution / belongs to 
H t (M),t<= R, if and only if 



11/11/ 




where Fourier coefficients Cj (/) of / are given by 

Cj(/) = (f,<Pj) = / fW 
Jm 

This L 2 — inner product can be also considered as a pairing between H_ t (M) 
and H t (M) and in this sense every element of i?_ t (M) can be identified with a 
continuous functional on H t (M). 

For a given finite family of pairwise different submanifolds {M v }i consider the 
following family of distributions 

(26) F V (J) = [ f 

which are well defined at least for functions in H e+d / 2 (M), e > 0. In particular, if 
M v = x v G M, then every F v is a Dirac measure 8 Xv v = 1, N, x v £ M. 
Note that distributions F v belong to H_ £ _ d / 2 {M) for any e > 0. 

Variational Problem. Given a sequence of complex numbers v = {tv}, 
v = 1,2, N, and a t > d/2 we consider the following variational problem: 

Find a function u from the space H t (M) which has the following properties: 

(1) F u {u) =v v ,v = 1, 2, N,v = {v„}, 

(2) u minimizes functional u — > ||(1 — £) t//2 tt||. 

We show that the solution to Variational problem exist and is unique for any 
t > to . Wc need the following Independence Assumption in order to determine the 
Fourier coefficients of the solution. 

Independence Assumption. There are functions £ C°°(M) such that 

(27) Fv{$n) = <W, 

where b~ vjl is the Kronecker delta. 

Note, that this assumption implies in particular that the functionals F v are 
linearly independent. Indeed, if we have that for certain coefficients 71, 72, 7at 

N 

I/=l 

then for any 1 < ji < N 

N 

The families of distributions that satisfy our condition include 
a) Finite families of 8 functionals and their derivatives. 
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b) Sets of integrals over submanifolds from a finite family of submanifolds of any 
codimension. 

The solution to the Variational Problem will be called a spline and will be de- 
noted as St(v). The set of all solutions for a fixed set of distributions F = {F v } and 
a fixed t will be denoted as S(F, t). 

Definition 4.2. Given a function f € H t (M) we will say that the unique spline s 
from S(F,t) interpolates f if 

F„(f)=F„(s). 
Such spline will be denoted as s t (/). 

From the point of view of the classical theory of variational splines it would be 
more natural to consider minimization of the functional 

\\£ t,2 u\\. 

However, in the case of a general compact manifolds it is easer to work with the op- 
erator 1 — C since this operator is invertible. The following existence and uniqueness 
theorem was proved in [34] . 

Theorem 4.3. The Variational Problem has a unique solution for any sequence 
of values (vi, v 2 , ...Ujv)- 

Proof. Consider the set 

P\KerF v = V t °(F) c H t (M),t > d/2, 

V 

of all functions in H t (M) such that for every 1 < v < N, F y (f) = 0. 

Given a sequence of complex numbers (i>i, v 2 , i>jv) the linear manifold 

V t (F,v u ...v N ),t>d/2 
of all functions / in H t {M) such that F„{f) = v Vl v = 1,...,N, is a shift of the 
closed subspace V t °(F), i.e. 

V t (F,v u ...,v N ) = V t °(F)+g, 

where g is any function from H t {M) such that F L ,(g) = v v , v = 1, 2, N. 

Consider the orthogonal projection go of g G H t (M) onto the space (F) with 
respect to the inner product in H t (M): 

(h,h) Ht(M) = ((i - cf'h, (i - £) t/2 / 2 ) L2(M) = 

/ (l-£)*/ 2 /i(l-£)*/2/ 2 . 

J M 

Note, that st(y) — g — go G Vt(F, vi, vn) is the unique solution of the Variational 
Problem. Indeed, to show that St(v) minimizes the functional 

on the set Vt(F, v±, ...,vn) we note that any function in Vt(F,Vi, ...,v jy) can be 
written in the form st(v) + h, where h 6 V t °(F). For such a function we have 



||(1 -£)*/ 2 ( St («) + /i)|| 2 = 
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||(1 - £) t/2 st(v)\\ 2 + 2 (s t (v), h) Ht{M) + ||(1 - C)^hf. 
Since s t (v) = g — go is orthogonal to V®(F) we obtain 

H(i - cy/ 2 ( St (v) + *h)\\ 2 = n(i - cy/ 2 St (v)\\ 2 + h 2 h(i - cy' 2 hf, h e v?{F), 

that shows that the function s t (v) is the minimizcr. □ 

The following criterion follows from the previous theorem [34] . 

Theorem 4.4. A function u £ H t (M) is a solution of the Variational Problem if 
and only if it is orthogonal to the subspace V®(F) and F v {u) = v v , v = 1,2, .... 

As a consequence we obtain that the set of all solutions of the Variational Prob- 
lem is linear. In particular, every spline st{v) £ S(F,t) has the following represen- 
tation through its values F v (st(v)) = v v ,v = 1, N, on X: 

N 



(28) **(«) = $>«A 



where F„(s t (v)) = v v , and I" £ S(F, t), v = 1,2, N, is so called Lagrangian spline 
that defined by conditions F^(l v ) = 5 vfl , fj, = 1, 2, N. 

The next Theorem gives the characteristic property of splines. 

Theorem 4.5. A function st(v) £ H t (M),t > d/2 is a solution of the Variational 
Problem if and only if it satisfies the following equation in the sense of distributions 

N 

(29) (l-Cy St (v) = Y / ^(st(v))'K. 

In other words, for any smooth ip 

r N 

((l-£y St (v)^) = / Q.-C) t 8t(vW = 52a v (8 t (v))F v (ffi). 

Proof. We already know that every solution of the Variational Problem is orthog- 
onal to V t °(F) in the Hilbert space H t (M) i.e. for any h £ V t °(F) 



(so) o = ( St (v),h) Ht(M) = {i-cy/ 2 s t {v){i-cy/*h. 

J M 

Note, that since we consider a finite family of pairwise different manifolds 
M v our Independence Assumption (|2T|) is satisfied: there arc functions d u £ 
C°°(M) such that 

(31) F„(0„) = S Vfl , 

where 6 vfi is the Kronecker delta. Thus, for any ip £ C yD (M) the function 

N 



belongs to V°(F) and because of j30 
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H t (M) 



In other words, 

JV 



JM „_! JM 



If we set 

a v (a t (v),4)= f [(1 -£)**(«)] ^> 

JM 

we obtain that (1 — £)*st(i>) is a distribution of the form 

N 

where 

So every solution of the variational problem is a solution of (|29[) . 

Note, that if C = {C^} is another Cfi°(M)— family for which F v ((fj_) = 5 vll , then 
we have the identity 

N 

Since distributions F v are linearly independent it implies that 

a^stiy),®) - a„(s t (v),() = 0. 

In other words, coefficients a v (st(v), t?) = a u (st(v),0 = a: u (st(v)) are independent 
of the choice of the family of functions 

Conversely, if u is a solution of then since .F„ belongs to the space H-t (M), 
and t > to > 0, the Regularity Theorem for elliptic operator (1 — Cf implies that 
u E H-t a+2t {M) C H t (M) and for any h e V t °(F) 

(«, fc) Bt(Jtf) = ((i - z^/v (i - cy' 2 h) = ((i - £)*«, h) 

N 

= Y^a v (u)F v (h) = Q, 

that shows that u is a the solution of the Variational Problem. □ 

The formula (|28jl represents spline through its values and Lagrangian splines. 
To obtain another representation of splines we will need what we call generalized 
fundamental solutions or generalized Green's functions £7*, which are solutions of 
the following distributional equations 

(32) (1 - CfEl = T v . 
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To find El we note that in the sense of distributions 

oo 

(33) F~ = J2FW<Pj 

3=0 

which shows that 

oo 

(34) Ei = ^(l + \ J )- t FJn)v j . 

3=0 

The following important fact holds. 

Theorem 4.6. Every spline s t (v) is a linear combination of the generalized fun- 
damental solutions 

N 

(35) s t (v) = Y / ^(st(v))El 

where scalar coefficients a v (st(v)) are the same as in Theorem \4-5\ 
Proof. According to Theorem 14.51 every spline St(v) is a solution of 

N 

(l-£) t « t (w) = 5^a v (s t (w))^. 
Thus, wc obtain the equality 

v 3 

that implies the desired representation. □ 

Note that so far we have used just the assumption that t > d/2. To get more 
information about s t (v) we will need a stronger assumption that t > d. 

Theorems 14.31 and 14.61 imply our main result concerning variational splines (see 

Theorem 4.7. If t > d, then for any given sequence of scalars v = {v u },is = 
1,2, ...N, the following statements are equivalent: 

1) st{v) is the solution to the Variational Problem; 

2) st{v) satisfies the equation \29\) in the sense of distributions 

N 

(36) (1 - Cy St (v) = J2<*Ast(v)jK,t > d, 
where u\(s t (v)\ ctpf(st{v)) form a solution of the N x N system 

N 

(37) ^ Pvn<Xv(st(v)) =Vn,fJ> = 1, —i N, 
and 

oo 

(38) ^ = ^(l + A J -)- t JU^)^fe'): &p J = -\<p j ; 

3=0 
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3) the Fourier series of St(v) has the following form 

oo 

(39) s t (v) = ^2cj(st(v))ipj, 
where 

N 

Remark 4.8. It is important to note that the system \31/\ ) is always solvable ac- 
cording to our uniqueness and existence result for the Variational Problem. 

Remark 4.9. It is also necessary to note that the series &38\) is absolutely con- 
vergent if t > d. Indeed, since functionals F v are continuous on the Sobolev space 
Hd/2+e{M) we obtain that for any normalized eig en function ipj which corresponds 
to the eigen value Xj the following inequality holds true 

\F V { V] )\ <C(M,F)\\(l-£) d /^\\ <C(M,F)(1 + A,y/ 4 , F = {F v }. 

So 

\Fj^)F tl { Vj )\ < C(M,F)(1 + A,) d / 2 , 

and 

|(1 + \ j T t F%M)F ll (<p j )\ < C(M,F)(1 + Aj)'* 0- *'. 
It is known [26] that the series 

E A 7 T < 

3 

which defines the £— function of an elliptic second order operator, converges if t > 
d/2. This implies absolute convergence of A38\) in the case t > d. 

One can show that splines provide an optimal approximations to sufficiently 
smooth functions. Namely let Q(F, f,t, K) be the set of all functions g in H t (M) 
such that 

(1) F„{g) = F v {f),v= 1,2,..., N, 

(2) \\g\\ t <K, for a real if > \\s t {f)\\ t . 

The set Q(F, f,t,K) is convex, bounded and closed. 

The following theorem (see |34j ) shows that splines provide an optimal approxi- 
mations to functions in Q(F, f, t, K). 

Theorem 4.10. The spline s t (f) is the symmetry center of Q(F, f,t,K). This 
means that for any g G Q(F, f, t, K) 

(40) \\s t (f) - g\\ t < X -diam Q(F, f, t, K). 

5. Approximate inversion of the group Radon transform using 
generalized variational interpolating splines 

We return to the Radon transform on a compact Lie group Q. In order to apply 
the general scheme developed in the previous section we select a finite number of 
pairs (x v , y v ) G Q x Q,v = 1,...,N, and introduce submanifolds j\A v = x v ,¥fy~ l C 
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Q. Note, that for any v = 1,...,N the dimension dim M. v = d v equals dim 3fc '. 
Next, for the set of scalars 



(41) v v =8$f{x v ,y v ) = f{x v hy,)^h= f(h)dh 

JJtf J M v 

we consider the following variational problem: for a t > ^dim Q find a function u 
in the space H t (G) which satisfies (|41[) and minimizes the functional 

u ^||(1-A 6 )*/ 2 k||. 



According to Theorem 14.71 the solution St = s t (v), v = {iv}, to this problem is 
given by the formula 

d„ 

(42) s t = J2 E 
where 

jv 

(43) cl } = cy(s t (u)) = (1 + A 2 )" 4 ^a^-^.j/,,), = a„(s t (v)). 

v=l 

Coefficients ai, ...,ajv are solutions of the following system 

(44) Pi^cti + ... + p Nfl a N = Vp, n = l,...,N. 

To determine matrix /? with entries /3 Vfjt , v, fi = 1, N, one only needs to find the 

data ^7r (g v ) and to compute quantities X]fj=i ^(^!j( x ^j y^))^(^ij( x M> The 
entries fi vyL arc given by the formulas 



(45) /V = E( 1 + A -)~* M ^ Xv ^ X ^ x ^ v ^' 

Tree M=i 

The appearing formulae for special applications will end up in well known calcu- 
lations of special function, since these arise naturally from representation theory. 
Implementing a fast algorithm for the solution of (|37ll , that is a standard problem, 
will give the solution of our variational problem for the Radon transform. 
Here we have 

(46) 

d„ d^ . . 

^2 @(mj(xi>,yv))@('Xij(xn,yii)) = E / ^iji^hyu 1 ) dh T7 l j(x f ,hy~ 1 ) dh 
i,j=l i,j=l Jjif Jjif 

(47) = y~] / / Kjiiyvhx' 1 ) dhTT^iy^hx^ 1 ) dh 

iJ=1 ■/ // Jjt 

(48) = trace(77^77(^)77 J ^77(a;^ 1 )7rj^7r(a; M )7rj^77(y A ; 1 )), 

hence we obtain function that is zonal in every component. A special case of that 
is the Addition theorem 13.31 for spherical harmonics. 
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In the case where JtiP is a massive subgroup of Q we find 

(49) W ( 

y v ))M{-K l:j {x ll ,y lJ )) = ir 11 (y„)ir 11 (x l/ 1 )ir 11 (x fl )ir 11 (y fl 1 ) 

i,j=l 

(50) =iru(y„x- 1 x ll y- 1 ), 
and hence 

(51) /V = E (! + X l r^niv^x^y- 1 ) 

where Q\ denotes the set of irreducible representations with rank tt^p = 1. 

Since 50(2) is a massive subgroup in 50(3) the last formula can be used in the 
case of the Radon transform on 50(3). 

Let {(xi, t/i), (xjv, J/jv)} be a set of pairs of points from Q. In what follows we 
have to assume that our Independence Assumption (|2"T)l holds. It takes now 
the following form: there are smooth functions <p\, </>_/v on Q with 

But it is obvious that for this condition to satisfy it is enough to assume that 
submanifolds M v = x^Jify" 1 C Q are pairwisc different (not necessarily disjoint). 
Let / be in a continuous function on (J, t > ^dimQ, v = {v v }^ where 

v v = f f- 

According to Definition 14.21 we use notation s t (f) = s t (v) for a function in H t (Q) 
such that for hA v = XyM'y^ 1 one has 

(52) {Mf) {x v ,y v ) = f f = [ s t (f) = (#*(/)) (x v ,y v ), 
and 

(53) l|si(/)l|i? t (S) ->-min. 

Theorem 5.1. Let {(x\, yx), (xn, yisr)} be a set of pairs of points from Q, such 
that submanifolds M. v = x v ^y~ x C Q, v = 1, N, are pairwise different. 

Given a continuous function f on Q and a t > ^dim Q the solution of I5£|) -(f? ^) 
is given by the formula CSD. The Fourier coefficients Ck(st(f)) of the solution are 
given by their matrix entries where a(st(f)) = (&v(st(f)))i S is the 

solution of {££]) with /3 £ R NxN given by [fify . 

The function St(/) S H t {Q) has the following properties: 

(1) s t (f) has the prescribed set of measurements 




*(/); 



(2) it minimizes the functional 

«-H|(l-A )*/ 2 u||; 
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(3) the solution is optimal in the sense that for every sufficiently large K > 
it is the symmetry center of the convex bounded closed set of all functions 
g in H t (G) with \\g\\t < K which have the same set of measurements 



/=/ 9- 

Let us turn to the case of 50(3). Since we have (l20l) 

47T 



(54) aT$(Z,T,) = T^TfM = ^—yi^yi^) 

it implies 

oo , . 2 2fe+l 

(55) ^ = ^1+^+1))-* (-^A EMw^w^w 

k=Q ^ ' i,j=l 

oo 

(56) = £(1 + k(k + l))-*ci(x v ■ y v )C*(xr ■ y J, 

k=0 

where we made use of Addition theorem 13.31 

At this point one uses standard methods to solve the last system, then to obtain 
coefficients ct v in (|4"4"]) . coefficients c^- in (|4"3")l and obtain the solution (|4"2"j) . 

Consider the constrained variational problem (|52"]) - ([5"3"| for Q = 50(3), J$? = 
50(2). Below we formulate a theorem which summarizes our results for 50(3). 

Theorem 5.2. Let {(xi,?/i), (x./v, 2Mr)} be a set of pairs of points from 50(3), 
such that submanifolds M. v = x lJ S0(2)y^ 1 C 50(3), v = 1,...,N, are pairwise 
different. 

For a continuous function f on Q, t > 3/2, and a vector (of measurements) 
v = ( v v)i where 

Vv = f, 

the solution of i52 \) - ^5^) is given by 

oo 2&+1 oo 

(57) St (f) = ^2Y, = Y, trace M s *(/)) Tfe ) . 

fc=Oi,j = l fc=0 

where T*. are the Wigner polynomials. The Fourier coefficients Ck(st(f)) of the 
solution are given by their matrix entries 

N 



(58) 4 (st(/) ) = (2fc+1)(1 4 ; fc(fc + i))* gM**(/)^)^), 



where a(st(f)) = (av(st(/)))f 6 M. N is the solution of 
(59) Pa(s t (f)) = /, 

with (3 £ M. NxN given by 



(60) = + k(k + l))-*cj(x v ■ y„)c| (*„ ■ y,), 

k=0 

l 

where Ci are the Gegenbauer polynomials. 

The function St(f) £ Ht(SO(3)) has the following properties: 
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(1) st(f) has the prescribed set of measurements v = (v v ) 1 at points {{x v ,y v ))i ; 

(2) it minimizes the functional 

u^\\(l-A so{3) ) t/2 ul 

(3) the solution |57j ) is optimal in the sense that for every sufficiently large K > 
it is the symmetry center of the convex bounded closed set of all functions 
g in H t (SO(3)) with \\g\\ t < K which have the same set of measurements 
v = (v u ) 1 at points {{x v ,y v ))i . 

6. Approximation by splines and a sampling theorem for Radon 
transform of bandlimited functions on 5*0(3) 

The Theorem 14.101 already shows that variational splines provide an optimal 
approximation tool. However, much stronger approximation theorems by splines 
hold in the case when the set of distributions F = {F„}^ is a set of delta functions 
on certain set of points {x v }^ of M (see [34] ). 

Definition 6.1. We will say that a finite set of points X p = (£i, £jv), p > 0, is 
a p-lattice, if 

1) The balls £?(£„, p/2) are disjoint. 

2) The balls £?(£„, p) form a cover of M of a fixed multiplicity Rm- 

The fact that the balls B(^ v ,p) form a cover of M of a fixed multiplicity Rm 
means that every point of M is covered by not more than Rm balls frm this fanily. 

The fact that for any manifold of bounded geometry there exist two constants 
Rm and po such that for any p < pa one can construct a p-lattice was proved in 

una, isa, [mi. 

We will need the following result for functions with " many" zeros from [50] , [52] , 

El, M- 

Lemma 6.2. There exist constants C(M) > 0,p(M) > such that for any p < 
p(M), any p-lattice X p = and for any f G H2d(M) such that f(£ v ) = for all 
£i, G X p , the following inequality holds true 

\\f\\<C(M)p 2d \\(l-C) d f\\,d = dimM. 

The next theorem extends the last estimate to higher Sobolev norms [53], [53]. 

Theorem 6.3. There exist constants C(M) > 0,p(M) > 0, such that for any 
< p < p(M), any p-lattice X p = {£„}, any smooth f which is zero on X p and 
any t > 

||(1 - CffW < {c{M) P 2d f m 11(1 - £f md+t f\\,t > 

for all m = 0, 1, .... 

It is important to stress, that the constant in the last inequality is growing 
exponentially when order of smoothness is increasing. 

Now we can formulate and prove our Approximation Theorem |33] , |34] . 

Theorem 6.4. There exist constants C(M) > 0, p(M) > such that for any 
< p < p(M), any p-lattice X p , any smooth function f and any t > the following 
inequality holds true 

||(1 - Cy(s 2md+t (f) -f)\\< {C{M)p 2 f md 11(1 - £) 2m+t /||, 
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for any m = 0, 1, .... 

Moreover, if t > d/2 + k then there exists a C(M,t) such that 

\\( S2 ^ + t(f)(x) - f(x))\\ ck{M) < (C(M,t)p 2 f md ||(1 - Cf md+t fl m = 0, 1, ... . 

As it was mentioned, using standard methods of interpolation theory and the 
fact that Bcsov spaces on manifolds are interpolation spaces (real interpolation 
method) between two Sobolev spaces the above inequalities can be extended to 
Besov norms. 

For a compact Lie group SO(3) and its closed subgroup SO (2) we consider 
S 2 x S 2 . According to Theorem O if / G H t (SO(S)) then &f G H^ +1/2 {S 2 x S 2 ). 
Also, integral of / over the manifold x u SO{2)y~ 1 is the value of Mf at {x Vl y y ). 
One can apply Theorem 14.71 to the manifold S 2 x S 2 and the set of distributions 
F = {F V ^} where 

F^ig) = gfr^vn), 9 eC(S 2 xS 2 ) 1 (x„, yfl ) e s 2 x s 2 

to construct corresponding intcrpolant s T (Mf) of £%f . Since dim (S 2 x S 2 ) = 4 for 
the intcrpolant s 2 ™+2 +t (3&f) £ H 2 m.+2 +t (S 2 x S 2 ) Approximation Theorem gives 
the estimate 

||(1 - 2A 52xS2 ) t ( S2m+2+t (^/) -^/)|| L2(s2xS2) < 
(C P 2 ) 2m+2 ||(1 - 2A s2xS2 ) 2m+2+t ^/|| L2(s2xS2) = 



(61) (C P 2 f \\Mf\\ H , 

for any m = 0, 1, .... 

Let s 2 ™+ 2+ i(^ 1 /) be orthogonal progection of s 2 ™+ 2+t (^?/) onto subspace 

Ranged = Hf /2 (S 2 x S 2 ). 

According to ((20)) it means that s 2 m + 2 +t(^'/) nas a representation of the form 

Applying (|20[) one more time we obtain that the following function defined on 5*0(3) 

s 2 ™+»+t(/)(0 = (W+t W)) (0 

has a representation 

s 2 - + » + t(/)(0 = EEir4W; m .^K). 

fc ij 

where are Wigner functions. 

Let us stress that functions S T (f) do not interpolate / in any sense. However, 
the following approximation results hold. 

Theorem 6.5. For t > 1/2 for sufficiently smooth functions one has the estimate 
||(1 - 4A so(3) ) t - 1 /4(S 2m+2+t(/) - /)|| La(SO(3)) < 

(62) {Cp 2 ) 2 ^ 2 ll#/llir am+a+t (s»xs») 
/or any m = 0, 1, .... 



x+2 



+2+t (S 2 xS 2 ), 



22 



SPLINES FOR RADON TRANSFORM ON COMPACT LIE GROUPS 



In particular, if a natural k satisfies the inequality t > k + 7/4, then 

(63) \\S 2 m+2 +t 

for any m = 0, 1, .... 

The proof of the first inequality follows from (59) and the fact that is 
continuous from H^{S 2 x S 2 ) to Sobolev i? t _i/ 2 (50(3)). The second inequality 
follows from the Sobolev Embedding Theorem. 

For an u > let us consider the space E u (SO(3)) of w-bandlimitcd functions on 
50(3) i.e. the span of all Wigncr functions with k(k + 1) < co. As the formulas 
(|17p and (j2"u) show the Radon transform of such function is w-bandlimited on 
S 2 x S 2 in the sense its Fourier expansion involves only functions y^y^ which are 
eigenfunctions of As^xS 2 with eigenvalue — k(k+l). Under our assumption about k 
the following Bernstein- type inequality holds for any function in the span of y^y^ 

||(1 - 2A s2xS2 y<%f\\ LHS2xS2) < (1 + 2w) T |W|| i2 ( S 2xS=) 
Note that if / G E w (S l O(3)) then the following functions are also bandlimited 

A;<u; ij 

s 2 ^ +t (f)(0 = E E ^^4W;™^-(£), 

fc<w ij 

Using our definition (|3.4[) of the norm in the space H t (S 2 x 5 2 ) we obtain the 
following refinement of Theorem 16.51 

Theorem 6.6. (Sampling Theorem For Radon Transform). 
For f S E w (<50(3)) and t > 1/4 one Aas ifte estimate 

||(1 - 4A SO(3) ) t - 1 /4(S 2m+2+t(/) - /)|U 2 (SO(3)) < 

(64) (Cp 2 (l + 2w)f m+2 (1 + 2iu) t ||^/|| 

for any m = 0, 1, .... 

in particular, if a natural k satisfies the inequality t > k + 7/4, then 

(65) ||S'2m+2 +t 

(/) - fWcHsom < {C P 2 (l + 2w)f m (1 + 2^*11^/11 

for any m = 0, 1, .... 

Let us remind that the function 5 , 2m+2 +t (/) was constructed by using only the 
values of the Radon transform Sif on a lattice of points on S 2 x S 2 . This Sampling 
Theorem shows that if 

1 

P < = 

y/UQ. + 2w) 

than one can have a complete reconstruction of w-bandlimitcd / when m goes to 
infinity by using only the values of its Radon transform &f on any fixed p-lattice 
of S 2 x S 2 . 
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